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We set up a zero recoil sum rule to constrain the form factors of the A;, —>■ transition. Our results 
are compared with the recent lattice calculation for these transitions. We find the same situation 
as in the case for B —>■ D*: The lattice results practically saturate the sum rules, leaving basically 
no room for excited states. 


I. INTRODUCTION 

The precise determination of the CKM matrix elements Vxt, x = c,u becomes increasingly important as an input 
for tests of the standard model at the precision level. Although lattice QCD as well as non-lattice methods - such 
as QCD sum rules - have made enormous progress, we are still facing a tension between determinations of Vxb from 
inclusive versus exclusive decays [T]. 

It is generally believed that \Vcb\ can currently be determined with the best precision via the inclusive decay 
B —>■ Xc^v [2 [3]. In this case one applies an operator product expansion (OPE) in terms of local operators, which 
sets up and expansion for the total rate, as well as for spectral moments, in powers of as and Agcn/mg, Q = b,c. 
This combined expansion seems to converge rapidly, giving us confidence in the precision of the method. 

On the other hand, exclusive decays also allow for a precise determination of I4h from the decays B —>■ D^*Hv by 
extrapolating to the point of maximal momentum transfer to the leptons At this point, heavy quark symmetries 
yield an absolute normalization of the form factors, and corrections to the form factor normalizations can be computed 
on the lattice HIS] as well as from QCD sum rules [SHni- 

The aforementioned tension between the inclusive and the exclusive determinations of |14f,| is driven by the lattice 
values for the form factor normalization for the B —> form factors. On the other hand, the anatomy of the b ^ c 

transition at zero recoil can be studied with zero-recoil sum rules, which hint at smaller values for the form factor 
normalizations and which are fully consistent with the inclusive determination. In particular, from the point of sum 
rules, the current lattice value would imply unexpectedly small contributions from the excited states [HlISj. 

More serious seems the problem with the determinations of \Vub\- The inclusive determination relies on a light-cone 
version of the OPE leading to the corresponding heavy mass expansion [T] . The hadronic input - the so called shape 
functions - are not well known (in particular at subleading order), and thus the resulting expansion leads to larger 
uncertainties compared to ones in the local OPE relevant for semileptonic b —b c decays. 

The exclusive determinations on Vub rely mainly on the channel B —> tt£P. For this decay, the form factors need 
to be computed either on the lattice m or estimated via light-cone sum rules mi. Using these form factors, which 
turn out to be consistent between the lattice and the QCD sum rules, a value of |I4b| can be extracted that is about 
three standard deviations smaller than the inclusive one. 

Since currently the exclusive determination of Vub rests mainly on a single channel, it is important to have an 
independent determination from an other channel. Since the purely leptonic decay B ^ £v suffers - even for the r 
lepton - from helicity suppression, the existing measurements oi B ^ tv are currently too imprecise to decide between 
the exclusive and inclusive value of Vub- This tension has also lead to sepeculations (see e.g. [12? ] that “new physics” 
is responsible for the effect, although right-handed currents have recently been excluded as an explanation |I8j . 

Recently the LHCb collaboration published a first measurement of the branching ratio of Af, —>■ p£v which is 

in principle precise enough to challenge determinations based on B ^ tt£v. However, this measurement is normalized 
to the branching ratio of Af, —J. Ac£v. Thus, the extraction of the ratio Y^ub/Vcb] requires the form factors to be 
calculated for both the Ab ^ p as well as for the Af, —J. Ac transition. This has been done recently on the lattice 
for both transitions with sufficient precision in m- Their results for the Af, — >■ p transitions compare favorably with 
light-cone sum rule calculations [16] . but the precision of these sum rules is intrinsically limited. 

In this work, we construct a zero-recoil sum rule (ZRSR) for the Af, —)■ Ac transitions, along the same lines as for 
the B —b D* form factor, see e.g. [8|. We shall investigate in this paper, if the tension present in the lattice calculation 
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versus the zero-recoil sum rule for the mesons persists for the case of the baryons. In the next section we formulate 
the zero-recoil sum rule for baryons and compute the necessary OPEs to the required level of precision. We apply this 
method to both the axial vector and the vector current, which eventually yields constraints for a subset of the from 
factors that describe the Af, —>■ Ac transitions. Finally we compare our results with the lattice values and conclude. 


II. ZERO RECOIL SUM RULE 

The sum rule at zero recoil ist set up in the same way as in the case for mesons [5] by considering the forward matrix 
element 


Trivq) = (Afc(P)|r{&(x)rc(x), c(0)r6(0)} |A,(P)) (1) 

where we shall discuss two possible choices of the currents T (g) T: 7 ^ (g (V x V) and 7^75 (g 7^75 [A x A). The 
normalization constants Ny are chosen to be Ny = 1 and Na = 3 for the two cases respectively. Furthermore, 
P = Ma^v is the momentum of the A^ baryon, from which we define the velocity v. 

We want to set up a sum rule at the kinematical point where the charm quark also moves with the same velocity 
V which is the point of zero-recoil transferred by the b ^ c transition. Thus we redefine the quark fields as 

b{x) = and c{x) = , ( 2 ) 

which suggests to define the parameter e = mt — rric — {v ■ q). We can then reparametrize the forward matrix element 
in terms of e, which leads to 

Trie) = ^J ("■")^(A,(P)|r{5„(x)rc„(x) c„( 0 )r 6 „( 0 )}|Afc(P)) (3) 

Since Ma,, — ~ mb — me-, the quantity e corresponds to the excitation energy of the intermediate charm states 

above the A^. The steps leading to the sum rule are formally as in [H], however, the relevant hadronic matrix elements 
will be different. Along the lines of [ 8 ] we define the contour integrals 

dn,r(£M) = j!> e"Tr(e)de, (4) 

\s\=Sm 

where the relevant contour is shown in figure 



Figure 1. Integration contour for the calculation of 7n,r(£M); the radius of the contour is em- Figure taken from [8]. 


Inserting a complete set of states, the lowest possible state is the A^ moving with velocity v, the higher states will 
excited states of the Ac but also non-resonant contributions such as AcTt or Dp, where the charmed hadron moves 
with velocity v. Looking first at the integral /o,r(£M) the lowest contribution thus is related to the square of the 
At, —> Ac matrix elements at zero recoil 


F 


1 


^(Ab(u, s)\by'yfj_Cy\Aciv, s')) (Ac(u, s')|c„ 7 '^&„|Af,(u, s)) 


(5) 


for the vector current, and 


G = ^ '^{Abiv, s)|6„7^75C«|Ac(r), s')) (Ac(u, s')\cy^''’^5by\Abiv, s)) 


(6) 
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Parameter 

mean value/ Icr interval 

unit 

prior 

source / comments 

quark-gluon coupling and quark masses 

as{mz) 

0.1184 ± 0.0007 

— 

gaussian @ 68% 

m 



nib {mb) 

4.18 ± 0.03 

GeV 

gaussian @ 68% 

m 



mc{mc) 

1.275 ± 0.025 

GeV 

gaussian @ 68% 

m 



hadronic matrix elements (nominal choice 


^2(1 GeV) 

0.50 ± 0.10 

GeV^ 

gaussian @ 68% 

see eq. ( 

24 

1 

pt(lGeV) 

0.17 ± 0.08 

GeV® 

gaussian @ 68% 

see eq. ( 

25 

1 


Table I. Summary of the prior PDFs used in the numeric analyses. The central values lead to = 0.75 GeV) = 4.62 GeV, 

= 0.75 GeV) = 1.20 GeV, and = 0.284. 


for the axial-vector current. 

We use the form factors for the > Ac transitions in the helicity basis, which is introduced in m- For the vector 
current they read 


{Ac{v\ s')\cjf_,b\Ab(v, s)) = MA,(w',s') 




i - {Ml - Mlf- 


S+ 

, 2Ma^Ma 

+fAw) ( 7^- {Vf, -h vl 




(7) 

( 8 ) 

(9) 


and for the axial vector current one has 


(Ac(u',s')|c757/i&|A6(t', s)) = -ua,('c',s')75 


go{w){MAt, + MaJ- 


+ 9+iw) 


Ma, - Ma, 


Ma,v^ + Ma^ - {Ml - Mlf- 


, , ,1 , 2Ma^Ma, , 

+5-l(w^) ( 7m +---(^M - ^m) 

'5+ 


uaAv,s) . 


( 10 ) 

( 11 ) 

( 12 ) 


In terms of the heavy hadron velocities v, v' and their scalar product w = vv' one finds q = Ma^v — Ma^ and 
g2 = Ml + Ml - 2Ma,Ma, w. In addition, we abbreviate 


s± = {AIa^ ± AIaJ^ - 


(13) 


With these definitions we obtain 


= l/o(w = 1)1^ , and 

G = ^ [2\gA{w = 1)|2 -h \g+{w = l)p] . 


(14) 

(15) 


The form factors f\ and g\, A = 0, -h, T, have been recently calculated on the lattice m, and are published in form 
of a handful of parameters, including their correlation matrix. Using their results for the form factors, the authors of 
FI obtain at the zero recoil point w = 1: 


F = 0.972 ±0.058, 


and G = 0.817 ± 0.044 . 


(16) 


In the rest of this paper we confront the above lattice results with the constraints obtained form the zero-recoil sum 
rule 
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The values shown here are taken from the arXiv version 3. 
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Figure 2. The result PDF (left) and the CDF (right) for the quantity Jo ,a as obtained from 10® random samples of the paameter 
space. We show the central 68% probability interval as the grey-shaded area. 


A. Axial Vector Sum Rule at Zero Recoil 


We start the discussion with the axial vector sum rule 

1 

Na 


Io,a{£m) = -^ X! s)| 5„7^75 Ci, |Xc{f)) (^c(t')l Cy^^j 5 by |At,(n, s)) 

Xc 1 

= G + GineliSAl) 


(17) 


In the above, G-meii^M) captures all inelastic contributions to the correlation function up to an energy em-, i-e., all 
contributions with excitation energies 0 < e < cm- Note that both terms G and CjneiC^M) are positive. We can 
therefore rewrite the sum rule as an upper bound for G: 


G < Io,a{£m) ■ 


(18) 


The left-hand side of eq. p7| ) can be evaluated in the OPE [5], and one obtains 

IoA^m) = - A/mAM,^i) - Af/„3(eM,Ai) +0(A;^ad/w^Ahad/m^) (19) 

where the perturbative contribution is the same as for the mesonic case 

= 0.75 GeV) = 0.970 ± 0.02, (20) 

which contains the as and the aj. corrections [5]. 

The power corrections differ from the mesonic results, since a priori the forward matrix elements for the Af, are 
different from the ones for the B mesons. Furthermore, for the A-like heavy baryons, the matrix elements of all 
the spin-triplet operators vanish. This is due to the fact that the light degrees of freedom do not have any angular 
momentum and thus cannot generate a chromomagnetic field. Hence, all matrix elements involving these operators - 
including AiA), p\sAb) ~ vanish. The non-perturbative power corrections for the baryonic case therefore read 


A 

A 


A 

A{A) 


2 ' 

l/m? 

4 

\ml ^ ml 

imbmc. 

A 

pUA) 

Pd{A) / 

' 1 3 

1 jmA 

Arn^ 

12mf, \ 

ml 


m{,mc 


( 21 ) 

( 22 ) 


The kinetic energy operator for the Af, baryon has been discussed in the context of the Af, baryon lifetime m- 
Using the spin-averaged heavy meson masses one obtains up to terms of order 1/m 


A{B)-A{A) 


2mf,mc 


C 


{{Mas 


MaJ - {Mb - Md)) (1 + 0{l/mA ■ 


(23) 


— m, 
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Figure 3. The result PDF (left) and CDF (right) of Ginei as obtained from 10® samples of both the OPE result for Iq^a and the 
lattice result for G. The grey-shaded area shows those results with Ginei < 0, which is unphysical. We find that the unphysical 
range of Ginei accumulates ~ 55% probability. 


The most recent results a of combined fit of the B-meson hadronic matrix eleme nts and Vch to the measured lepton- 


energy moments in B —)• yield li^(B) = (0.47 ± 0.07) GeV 


Using eq. (241 this translates to 


Ai^(Ab) = (0.50 ±0.10) GeV^ 


(24) 


where we increase the uncertainty to account for the lack of 1/m^ terms. 

Given the small difference between the kinetic energy parameters of baryons and mesons, we use also for the Darwin 
term of the Af, the same value as for the B-meson. The mesonic matrix element is obtained in [20) : for the Af, we use 


e same central value and increase 

the uncertainty by a factor of two. 



p%Ab) - (0.17 ± 0.08) GeV® . 

(25) 

Using these numbers, we obtain 




{eM=P = 0.75 GeV) = 0.108, 

(26) 


(em = At = 0.75 GeV) = 0.028. 

(27) 


We note that is about 20% larger than for the mesonic case, while Ai/,^3 for the A^ baryon yields numerically 

the same result as for the meson. The above results shall only be illustrative, and have been obtained for our default 
choice of input parameter. 


For a more thorough numerical study, we use and extend EOS This allows us to carry out a Bayesian 

uncertainty propagation based on Monte Garlo techniques. We choose prior probability density functions (PDFs) 
for all input parameters based on the principle of maximum entropy [22) . We use Gaussian distributions throughout 
this work, since in all cases the mean and variance of the parameters are known. For a summary of the PDFs, see 
table 1^ We draw 10® random samples from B(/o,a), the PDF of our quantitiy of interest. The result PDF and the 
corresponding Cumulative Probability Density Function (CDF) are shown in figure]^ For our choice of the prior 
PDFs, the result is a gaussian distribution to very good accuracy, with skewness —0.08 and excess kurtosis of —0.04. 
From the result PDF we obtain the mode and the central 68% probability interval 

IoA^m = f^ = 0.75) = 0.8111°;°26 • (28) 

Following eq. (|17|), we can also compute the inelastic contributions Ginei from our nominal results for Iq^a and the 
lattice results foiM7. Using 10® samples of both quantities, we obtain the PDF and CDF for the quantity Ginei as 
shown in figure Again, the PDF is approximately gaussian with skewness of about —O.Ol and excess kurtosis of 
about —0.006. We obtain the mode of the distribution and the central 68% probability interval as 

Ginei = -0.005^°;O« . (29) 


Roughly 55% of the samples of Ginei turn out to be unphysical, since they are negative. Thus we conclude from this 
statistical analysis that the situation for the A;, —)■ Ac is very similar as for the B —)■ D* case: The lattice results for 
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Figure 4. The result PDF (left) and CDF (right) for the quantity Io,v as obtained from 10® random samples of the paameter 
space. We show the central 68% probability interval as the grey-shaded area. 


the form factors at zero recoil saturate the corresponding zero-recoil sum rule by a very large degree, leaving almost 
no room for inelastic contributions. In fact, compared to the mesonic case, the situation seems to be even worse, since 
the central value obtained from the lattice eq. (16 1 exceeds the central value for our upper bound. Furthermore, for 
the mesonic case, one may estimate the inelastic contributions, which turn out to be sizable. This in turn implies 
that the zero-recoil sum rules would predict a smaller value for the form factors. Unfortunately, the estimates in the 
mesonic case rely on the so-called BPS limit, which cannot be used in the case of baryons. Since an estimate of the 
inelastic contributions in the case of the A;, requires (possible even model dependent) input, we will not discuss this 
in the present paper. 


B. Vector Sum Rule at Zero Recoil 


The vector sum rule is obtained from eq. (|^ by inserting T 0 T = 7 ^ (g) and Ny = 1, 

1 

Nv 


Io,v{£m) = ^ {Ab{v,s)\by^^Cy\Xc{v)) {Xc{v)\cyj>"by\Ab{v,s)) 

Xci 

= F + Finel(eM) ■ 


(30) 


Analogous to the axial vector current, Finei(£M) captures all inelastic contributions to the correlation function with 
excitation energies less than em, i-S., all contributions with excitation energies 0 < e < em- Again, F and Finei(£M) 
are positive, and we can therefore rewrite the sum rule as an upper bound for the term F : 


F < loyisM) ■ 


The OPE result for the left-hand side of eq. (301 reads 


Io,v{£m) — (sM , fi) — AU„2(eM,Ai) — p) + C>(Ahad/^&i Ahad/w^) 


(31) 


(32) 


where the perturbative contribution has been evaluated to order as in [23] . For the central values of the input 
parameters we obtain 


ev"\£M = M = 0.75 GeV) = 1.03t°;°? , 

where the uncertainty is estimated form a variation of the scale 0 < /r < 1.5 GeV. 
The nonperturbative corrections have been given in [6l[23] 


= 


lil[Ab) ( 1 1 


A 


l/r; 




rUc rUb 

1 1 

rric mb 


1 1 
me mb 


(33) 

(34) 

(35) 
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Figure 5. The PDF and the CDF of Flnei as obtained from 10® samples of both the OPE result for /o,v and the lattice result 
for F. The grey-shaded area shows those results with Flnei < 0, which is unphysical. We find that the unphysical range of Flnei 
accumulates ~ 55% probability. 


and reflect the fact that the vector current is conserved in the limit mb = me- 
Inserting the central values for the hadronic matrix elements, we obtain 

(36) 

(37) 


{eM = tJ- = 0.75 GeV) = 0.047, 
(em = At = 0.75 GeV) = 0.017. 


As before, these results are only meant as an illustration, and we repeat the statistical procedure as outlined in 
We obtain for the mode and central 68% probability interval of the result PDF for loy 


section 


IIA 


Ioy{eM = At = 0.75) = 0.965 ± 0.013, (38) 

based on 10® samples. We display the resulting PDF and CDF for Igy in figure We compute the inelastic 
contribution as well - just as before in the case of the axialvector current - and obtain 


Finei = -0.010t°;°®i. (39) 

as the mode and central uncertainty interval at 68% probability; see figure for the respective result PDF and CDF. 
We further find that ^ 55% of the drawn samples are unphysical, i.e., they show a negative inelastic contribution. 

Thus our findings are qualitatively the same as in the case of the axial current: The lattice result for the scalar 
vector form factor /o at the non-recoil point again saturates the the zero-recoil sum rule to a very large degree, leaving 
also for this case almost no room for inelastic contributions. 


III. DISCUSSION AND CONCLUSION 

The determination of CKM matrix elements from exclusive semileptonic decays requires reliable calculations for the 
form factors describing the corresponding hadronic transition. Since the form factors are genuinely non-perturbative, 
the only known “ab initio” calculational method is lattice QGD. The progress in this held made in the last years in 
the construction of efheient algorithms as well as the increasing computing power has turned lattice calculation of 
form factors into an indispensable tool in flavor physics. 

However, despite this progress it is important to perform checks of the lattice results from “continuum” methods. 
One of these methods are QGD sum rules. On the one hand they are hrmly rooted in QGD, on the other hand they 
allow for a detailed study of the “anatomy” of the results obtained e.g. for form factors. It has to be clear that a 
QGD sum rule can never make a precision prediction for a hadronic quantity, since the method is intrinsically limited 
to a level of a few ten percent. 

Nevertheless, QGD sum rules can serve to validate results obtained from other methods, e.g from lattice QGD. In 
particular, the zero-recoil sum rules can give a hint on the sizes of the from factors at the non-recoil point; in case of 
the B —> D* transition one can combine the zero-recoil sum rule with an estimate for the inelastic contributions to 
actually estimate the form factor itself. 










In the analysis presented here we have shown that the lattice results m for the Af, —)■ Ac transition form factors 
saturate the zero-recoil sum rule to a large extent. In fact, we found that the central values for the lattice results 
exceed the sum rule’s upper bounds, leaving practically no room for any inelastic contribution. This seems to be the 
case for both the axial-vector as well as for the vector current. 

In fact, the degree of saturation of the sum rule for the A;, —>■ Ac seems to be higher than for the B D* transition, 
where the lattice value for the form factor at zero recoil still leaves room for a (too?) small inelastic contribution. 
Unfortunately, the inelastic contributions for the baryonic case are harder to estimate than in the mesonic case; any 
estimate of the inelastic contributions for the baryons would require (probably model-dependent) additional input. 
We leave the discussion of this to future work. 
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